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Abstract 

We study radial perturbations of general relativistic stars with elastic matter sources. We 
find that these perturbations are governed by a second order differential equation which, along 
with the boundary conditions, defines a Sturm-Liouville type problem that determines the eigen- 
frequencies. Although some complications arise compared to the perfect fiuid case, leading us to 
consider a generalisation of the standard form of the Sturm-Liouville equation, the main results 
of Sturm-Liouville theory remain unaltered. As an important consequence we conclude that the 
mass-radius curve for a one-parameter sequence of regular equilibrium models belonging to some 
particular equation of state can be used in the same well-known way as in the perfect fluid case, 
at least if the energy density and the tangential pressure of the background solutions are contin- 
uous. In particular we find that the fundamental mode frequency has a zero for the maximum 
mass stars of the models with solid crusts considered in Paper I of this series. 

PACS: 04.40.Dg, 97.10.Sj, 97.60.Jd 



1 Introduction 

This is the second paper in a series that intends to bring the relativistic theory of clastic matter sources 
to a point where it can be conveniently used in numerical studies of compact stars. In the first paperjT], 
hereinafter referred to as Paper I, we reviewed and extended the theory as well as applied it to static 
spherically symmetric (SSS) configurations. As an illustration, we constructed stellar models with elastic 
crusts or cores numerically using a relativistic polytrope as the unsheared part of the equation of state, while 
the shear modulus was assumed to be proportional to the pressure. 

In this paper we study adiabatic radial perturbations of general relativistic SSS stars with elastic matter 
sources. By an adiabatic perturbation we effectively mean that the cold matter equation of state that is 
used to model the eqillibrium configuration is assumed to be valid also for the perturbed, dynamic, system. 
Thus we limit our analysis to situations where the time scale of the perturbations is much smaller than any 
time scale associated with reactions within the medium. It should be mentioned here that this need not 
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be a good approximation for the matter in the deep interior of neutron stars, for a discussion see [2] and 
references therein. 

The study of steUar oscillations dates back to 1918 when Eddington wrote down the equation for adiabatic 
radial oscillations of Newtonian perfect fluid stars[3|. The first landmarks in the study of the corresponding 
general relativistic case were the works of Chandrasekhar^E] including his famous result on the dynamical 
instability of white dwarfs|Sj. Soon thereafter Tooper[7] wrote a very thorough paper on the equilibrium 
and stabillity of stars obeying a relativistic polytropic equation of state. One also has to mention the 
book by Harrison et al. |H]. Since then, a huge pile of papers have been written on the subject but 
we make no attempt to review the historical progress of the field here. Some more recent references to 
previous works can be found in the nice paper by Kokkotas and Ruoffjni where a large set of frequencies is 
calculated for realistic tabulated equations of state as well as classical polytropes. Radial perturbations of 
SSS models with anisotropic pressures have been studied before by Hillebrandt and Steinmetz^UI as well 
as Dev and Gleiser[TT]. However, in those works ad hoc ansatze are made for how the matter responds to 
the perturbations. Their results give qualitative insights into how pressure isotropy affects the stability of 
SSS models, but when using a particular matter description one wants to let that description govern the 
equations. The Einstein equations for dynamical spherically symmetric elastic models was notably studied 
by Magli in jl2) . but that work was not directly concerned with small oscillations around static models. 

This paper is organised as follows; In section [5] we consider radial perturbations of generic spherical 
matter sources using a coordinate independent formalism. The only restriction we put in at this time is 
that the stress-energy tensor should be diagonalizable in an orthonormal frame. This is obviously a very 
mild restriction, since most normal matter sources are believed to be of this type. We take the somewhat 
non-standard approach of using the Lagrangian (or comoving) gauge from the outset rather than starting out 
from the Euler (or Schwarzschild) gauge only to transform to the Lagrangian one along the wavjl3j. We then 
specialise to a static background which of course simplifies the problem considerably. Taking as our strategy 
to eliminate all occurrences of the perturbation of the two-metric orthogonal to the SO (3) orbits we end up 
with two equations relating the perturbations of the matter degrees of freedom to the perturbation of the 
Schwarzschild radius. These equations, being derived without specifying any particular matter type, should 
prove useful for just about any situation where radial perturbations of an SSS background are considered. 
Indeed, every SSS metric corresponds to a stress-energy tensor which is diagonal in an orthonormal frame 
and completely described by the energy density, the radial pressure and the tangential pressure as measured 
by a static observer. When radially perturbed, the stress-energy tensor in fact always stays diagonalizable 
(to first order) with respect to an orthonormal frame unless the energy density and the radial pressure sum 
to zero, in which case the diagonalizability may or may not continue to hold. For any physically reasonable 
matter type that is used to model the interiors of stars, this implies no restriction. 

In section |21 we consider the boundary conditions using the general Israel method of matching the first 
and second fundamental forms. For the case of a static background on which the energy density and the 
tangential pressure are continuous (the background radial pressure must always be continuous) we find that 
these conditions boil down to the matching of the Lagrangian perturbations of the Schwarzschild radius and 
the radial pressure. If, however, first order phase transitions in the energy density or tangential pressure are 
present the matching is not uniquely determined by the Israel junction conditions. Thus, in such situations, 
one needs to make physical assumptions about the nature of the phase transitions. At the surface of the star 
we find that the Lagrangian perturbation of the radial pressure should vanish, whether or not continuity of 
the energy density and the tangential pressure hold there. 

In section 0] we finally specialise to elastic matter using the formulation presented in Paper I. We find 
that, by giving an elastic equation of state as an input, the problem reduces to a single second order ordinary 
differential equation much in the same way as for perfect fluids. Unlike the perfect fluid case, however, the 
problem is not of the obvious Sturm-Liouville type, which is closely related to the fact that the Lagrangian 
perturbation of the radial pressure can not be written as a background function times the first derivative 
of the perturbation function that for perfect fiuids is the unknown variable in the Sturm-Liouville equation. 
The obvious way of proving that one has a Sturm-Liouville problem at hand is to make a transformation 
to the standard form of the Sturm-Liouville equation. This can indeed be very easily done, but when doing 
so one introduces the derivative of a background function which may have jump discontinuites. For this 
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reason we were not absolutely sure that our perturbation equation with the imposed boundary conditions 
really defines a system with all the standard Sturm-Liouville properties, so we go through some lengths 
to prove that so is in fact the case under quite general assumptions. Some parts of the proof concerning 
the qualitative properties of the Priifer type phase function of the system are somewhat technical and are 
therefore put in an appendix at the end of the paper. We conclude this section with a discussion on how the 
static approach to stability generalises to the case of elastic matter. 

Regarding the notation used in the present paper, we stick to the notation of paper I as far as possible. In 
addition, we use A and S to denote Lagrangian and Eulerian perturbation operators, respectively. Whenever 
one of these operators appears in an equation, it is to be understood that the equation is evaluated on the 
background spacetime. 

2 Radial perturbations of static spherically symmetric spacetimes 

Our approach to radial perturbations is to start out by writing the full nonlinear Einstein equations for 
a general spherically symmetric system in covariant two-dimensional language. To this end we denote the 
metric on the SO (3) orbits by tab {t for tangential) and the Lorentzian 2- metric on the orthogonal spaces 
M± by jab- Furthermore, the curvature radius of tab will be denoted by r while the connection, Gaussian 
curvature and volume form of jab will be referred to as T>a, K and Cabi respectively. We also introduce the 
mass function m defined by 

1- — =I?VI?,r. (1) 

r 

The four-dimensional Einstein tensor Gab can now be written as (c/. appendix A of ref. |14| ) 



^ Ga. 



b ■- 



- - j jab\ + {V'T^cT - rK)tab. (2) 
Thus the stress-energy tensor satisfying Einstein's equations Gab = uTab must be of the form 

Tab ~ %ib + Pttab, ^abt^ c = ^, (3) 

where pt will be referred to as the tangential pressure. The Einstein equations obviously split into 

-VaVbT + [v^V.r - ^) jab = \nr Tab (4) 
VVaT - rK = nrpt (5) 
Taking the trace of equation and substituting the result back we find 

VaVbr = ^jab - \i^r{Tab - Tjab) (6) 



r 



9m 

K = — + \K{T-2pt) (7) 



It is worth noting that a contraction of eq. © with D^r gives the following expression for the gradient of 
the mass function: 

Vam = ^Kr^iTa'' ~ Tja'')Vbr (8) 
The matter equations of motion VbT''a — can be written in two-dimensional language as 

r-^Vbir^T^,) - 2r-'ptVar = VbT\ + 2r-\Ta'' - ptja'')Vbr = 0. (9) 

The fact that K is supposed to be the Gaussian curvature associated with the connection Va can be expressed 
as 

e''''VaVbVc^ Kec^Va, (10) 
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for any one- form Va- Taking Va — "DaV and using eq. © we may rewrite the left hand side of eq. ()lUfl as 

e'^^VaV^V.r = + W\ e,-Var - \nr-^Vk{r^T\)e,\ (11) 



Setting this expression equal to Kec"''Dar we find that eq. is in fact implied by eqs. © and ©• We 
shall use this fact to let the the latter two equations provide the full Einstein equations for the spherically 
symmetric system. 

We shall now assume that Tab is diagonalisable in an orthonormal frame, viz 

jab = -UaUb + ran (12) 
Tab^ pUaUb+PrTaTb, (13) 

where —u°'Ua = r°'ra = 1, u^Ta — 0. We refer to this frame as the comoving frame and the scalars p and 
Pr as the energy density and the radial pressure, respectively. The symmetric connection Va is completely 
determined by the single connection one-form 

r, = r^'VaUb = -v^Van. (14) 

In order to simplify the notation we introduce timelike and spacelike derivatives according to 

A u'^VaA (15) 
A' r''VaA (16) 

for any scalar or tensor A. The connection one-form is in turn completely determined by the accelerations 
of the two frame vectors. Indeed, it can be expressed as 

Ta = -{r^Ub)ua - {u^r'^)ra, (17) 

where we note that the accelerations can be interpreted either as two-dimensional or four-dimensional objects 
since iib = u'^VaUb — u'^WaUb and = r'^Varb = r°'Warb- Moreover we denote the frame projections of the 
connection one-form by 

To u'^Ta = r'^Ua, (18) 
Ti := r'^Ta = -u'^C (19) 

With the dot-prime notation the definition of the mass function |^ takes the form 

_^2+/2^1_^. (20) 

r 

Projecting eq. © with r°'r'' and u'^^r''^ we find 

Tor' = - ^nrpr (21) 



r" - Tif = ^ - ^Krp (22) 



r(''^-^(ror + rir')=0 (23) 

where the symmetrization of the mixed derivatives of r, r^''\ is needed because the derivatives do not 
commute. The projections of the matter equations ^ take the forms, 

p+ip + Pr)Ti + 2ip + pt)- =0, (24) 

r 

p', + {p + Pr)To + 2{pr - pt)- ^ 0. (25) 

r 
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These last five scalar equations comprise the full Einstein equations in spherical symmetry. For completeness 
we also write down the projections of the mass gradient equation (jSj). They are simply 

m = —^Kr^prV (26) 
m' = ^Kr^pr'. (27) 

If the spacetime is static it is easy to see that eqs. (|21|) . (|22|l and (|25|1 . when expressed in Schwarzschild 
coordinates, reduce to the well-known Einstein equations for a general SSS spacetime 15 , with eq. (|25|l then 
being the generalised TOV equation of hydrostatic equilibrium when anisotropic stresses are allowed for. 
Furthermore eqs. (|^ . and reduce to identities in the static case. 

2.1 Perturbation formalism 

We shall apply a comoving, or Lagrangian gauge, in which AA can be interpreted as the perturbation of a 
quantity A as measured by an observer moving along with the four- velocity of the matter. This is compatible 
with 

AUa=l{Ua, (28) 

Ar, =7^r„ (29) 

for some scalars U and TZ. We take this to be the definition of the gauge. ^ In contrast, one could also 
make the Schwarzschild, or Eulerian, gauge choice. Then 6A denotes the perturbation of A as measured 
with respect to an observer sitting at a fixed point in the Schwarzschild coordinate grid. To discuss the 
connection between the two gauge choices it is convenient to introduce the Schwarzschild frame given by 

jab = -UaUb + rah, (30) 

where 

u^Var^O, = f-^f, = 1, u"fa = 0. (31) 

The relation between the comoving and the Schwarzschild frame is uniquely specified by the boost velocity 
V that relates the two frames according to 

= j{u'' - vr") (32) 
= jir" - vu"). (33) 

where 7 = (1 — w^)~^/^. The Schwarzschild gauge can be defined by 

6r^0, dua^Uua (34) 

which implies 5fa = T'^fa-, for some scalars U and TZ. It is common practise to start out from the Schwarzschild 
gauge and subsequently transform to the comoving gauge. We shall instead use the comoving gauge from the 
outset and discuss the connection to the Schwarzschild gauge afterwards. Now, let Cab be the background 
tensor which can be viewed as the Lagrangian perturbation of the two-dimensional Christoffel symbols in 
the sense that, for all one- forms Va, 

A{VaVb) = VaiAvb) - VcC'ab- (35) 

One can show that this tensor satisfies 

r^cC^a = VaTl (36) 

u''u,C%a = -VaU, (37) 



^strictly speaking this is not a proper comoving gauge which would instead be defined by Au" = 0. 
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from which it follows that 



Perturbing eq. (pn|l we find 



AFo = w'-ro7e, 

AFi = n- TiU. 



while perturbing eqs. (|?T|l - (|^ yields 

i - ToC -fU-r'U' - 2{rU - Tar'TZ) 



-A - + — ^ - InirApr + PrO 
r \ r / 



C - Tie - rTt - r'W - 2{r"n - T^rU) ^Ia(-) ^K{rAp + p^) 

r \ r / r 



(38) 
(39) 

(40) 

(41) 
(42) 
(43) 



(Ap) • + Ti ( Ap + Ap.) + (p + p,.)^ + y ( Ap + Apt) + 2(p + ft ) ""^ ^ (44) 
(Ap,)' + ro(Ap + Aft) + (p + ft.)Z^' + ^(Aft - Aft) + 2(ft - ft) '^^ ~/ ^ = (45) 

This is as far as we intend to proceed in full generality. In the next subsection we specialise to the case of a 
static background with the application to radially oscillating stars in mind. However it would be interesting 
to use the above set of equations to study perturbations of collapsing or exploding configurations. 

2.2 Static background 

With a static background the equations given above simplify considerably. The unperturbed four-velocity 
w° will then by necessity (unless p + pr = 0) be aligned with the static Killing vector t", i.e. 



(46) 



where v can be interpreted as the gravitational potential. This means of course that 5 = for all unperturbed 
scalars S. It also follows that Fq is simply given by 



To = ^\ 



(47) 



whereas Fi vanishes identically since the acceleration r'^ does. Consider first equation H43(l which simplifies 
to 
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ryl 



(48) 



Hereinafter we shall set all static integration functions to zero since they are only interesting when perturbing 
a static solution to a static neighbour. Hence eq. H48() integrates to 



Likewise, eq. (I44|) is also a total dot derivative whose integrated version is 



Ap+ -:{p + Pr 



r 



where q = ^{pt ~ Pr) and equation (|48|l has been used. We also solve eq. 141|) for W', with the result 



W = -\\KrApr+'i+v'^' 



1 I i I I \ III' 

Upr - J.' 2i.' + - - - 



(49) 



(50) 



(51) 
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Using eqs. (|49ll . (|50|l as well as the background equations we may rewrite (|45ll according to 



-e " r \ — e Ap, 
r' 



P + Pr 



6 

r 



4r' 



(52) 



r'Aq } = 0. 



These last four equations govern the dynamic radial perturbations of any SSS background and regardless 
of the matter sources, provided that the stress-energy tensor is diagonalizable in an orthonormal frame. 
As already mentioned in the introduction, the condition of diagonalizability is very unrestrictive since it 
is automatically fulfilled if p + 7^ holds on the background. For many matter types the equation of 
state only imposes relations between the perturbed matter quantities Ap, Ap^, Aq as well as the perturbed 
Schwarzschild radius ^. In such cases the two equations (|5n| and (|5^ . from which all reference to the 
perturbation of the two-metric jah has been eliminated, will give all the information needed. 

Let us now briefly discuss how to transform a perturbed quantity from Lagrangian to Eulerian gauge, in 
case one would like to do such a transformation. If the quantity is a tensor field the relation is given by a 
vector field X°' according to 

(5 = A-£x- (53) 

Since we are restricting ourselves to a static background the Lie derivative will only be applied to static 
tensor fields. In fact we only need to apply it to static scalars which means that it is sufficient to determine 
the component = raX'^ which according to the gauge condition (5r = is simply found from 



= 5r = e - Cxr = C - X^r' 



(54) 



So, for practical purposes we may set = {^/r') r". Consequently the perturbation operators, when acting 
on scalars, are related according to the well-known formula 



, 9 

d = A-^—, 
or 



(55) 



if Schwarzschild coordinates are used. 



3 Boundary conditions 

Quite generally, we want to consider the junction conditions across a hypersurface H of the perturbed 
spacetime which is given by 

F = constant, (56) 

where is a spherically symmetric function. We assume that the surface is timelike so that it has a spacelike 
unit normal which can be written as 

n-a = Jnira ~ VHUa), 7h = (1 - uj)^^/^, (57) 

where vh — F/F' can be interpreted as a boost velocity. This boost relates the two orthonormal frames 
{u°',r'^} and {w°',n'^}, where w°' = 7ij(u" — VHr""). The Israel junction conditions imply that the first and 
second fundamental forms, qab — 9ab — naUb and Kab = qa'^qb'^^cnd, should be matched continuously across 
H. These quantities may be expressed as 

Qab = -WaWb + tab (58) 

Kab = iw''w'^Vcnd)waWb - {r-^n''Vcr)tab, (59) 

where we have used the relation 

S/aUb = VaUb - {r^^n'"Dcr)tab- (60) 
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This directly implies continuity of Wa as well as of the scalar quantities r, n°''Daf and w'^w'^Vcnd- It will 
prove convenient to replace the last of these quantities by 



where Rh is the Ricci scalar of H. The continuity of J follows from the relation (c/. jl^) 

Rh = -2Gabn''n'' + {K\f - X^^'i^afc- 

Using now that Gatn'^n'' and Rh can be expressed as 

Gabn'^n'' = KJniPr + Vh p) 

Rh = ir-^w'^Vaiw^Vbr) + 2r-^[{w''Varf + 1], 

we find that the Lagrangian perturbations of Wa, r, n°''Daf and J are 

(AF)' 

AWq ^UUa-\ ra 

Ar = ^ 

A(n"2?ar) = - r'TZ = v'S, 



A J = K,/S.pr 



2r- 



£ — r 



(61) 

(62) 

(63) 
(64) 



(65) 

(66) 
(67) 

(68) 



where we have directly specialised to the case of a static background, which for practical purposes implies 
that H to zeroth order coincides with a surface r — constant and that vh vanishes when unperturbed. To 
obtain perturbed scalars that are to be matched continuously over H we make a gauge transformation to 
a perturbation operator Sh for which ShF ~ 0, thus corresponding to an observer that stays on H . This 
gauge transformation is given by 



AF 



(69) 



where has no a priori restriction. Since the unperturbed values of Wa, r, n'^Var and J are Uo, r, r' and 
npr, it now follows that 



SnWa 



u-.'^-ix^) 



F' 



(X" 



^Hr = C-—r 

AF AF 

Snin'^'Dar) = i^'£_ - r" = i^'Snr + \nr{p + pr) 

Oh J = KApr + 2r (dHr) - np^ = n y Apr ~ —r j 



F' 



v'5H{n''Var)-{u')HHr+{5Hr)- 



(70) 

(71) 
(72) 

(73) 



where background equations have been used. Regardless of the function F, we can clearly make r^-dnWa 
continuous by choosing the function X^ appropriately. We can also make u'^SnWa continuous by making 
an appropriate matching of U across the boundary, which is possible since the perturbation equations only 
determines W by eq. (|5T|l . Thus the continuity of SnWa can be assumed to have been taken care of. We now 
note that the background quantities r, r', pr, v and v' are required to be continuous, where the continuity 
of Pr and v' are equivalent. Moreover, when assuming a time dependence of the type e^*"* for all perturbed 
scalars S, we can also use that continuity for S implies continuity for dot derivatives of 5 (unless w = 0). It 
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thus follows that the junction conditions for the radially perturbed solution imply continuity for the three 
quantities 



S2 = {p+Vr)T 



(74) 
(75) 




(76) 



r 



where 



Af 
F'/r'' 



(77) 



These junction conditions in general depend on the choice of the function F that defines the boundary hy- 
persurface H . Some important special cases are worth mentioning here. If p is continuous on the background 
it follows from the continuity of S2 that T is continuous which in turn implies that ^ is. If p as well as the 
pressure anisotropy scalar q are continuous on the background, then so should both ^ and Ap^ be. Moreover 
Apr is always continuous over a phase boundary if q vanishes on both side of that boundary. This is of 
course always the case for perfect fluids. 

For perfect fluids it interesting to consider the situation when the background solution is such that p 
has a jump discontinuity occuring at some pressure pi. In this case the Israel junction conditions do not 
determine the quantity T across the boundary. Therefore it must instead be determined by the microphysics 
that governs the state of matter. Then there are two especially simple cases where the oscillation modes can 
be computed without a detailed knowledge of that microphysics. 

The first alternative is to assume that the reaction rate is much faster than the timescale associated with 
the perturbation. Then matter can always be considered as being in thermodynamic equilibrium which mean 
that the phase transition can be viewed as occuring at pi in the perturbed star also. One should thence take 
the function F to be the pressure p Vr\ which gives 



Inserting this into eqs. (|74|I - H76(I leads to continuity of two quantities only, namely f — [p' /r')^^lS.p and Ap. 
In this case particles that are close to the phase boundary are allowed to cross that boundary, meaning that 
observers that are comoving with these particles will experience that the energy density changes discontinu- 
ously. As a consistency check of these boundary conditions, we have numerically studied radial oscillations 
of models with a one-parameter family of perfect fluid equations of state p = p{p, a) for which p develops 
a jump discontinuity as a function of p when the parameter a tends to a certain value ai. We found that 
whereas ^ also develops a jump discontinuity as a — * ai, neither ^ — (p' /r')^^Ap nor Ap do. 

The second alternative is to set AF — 0, which implies T = 0. This has the interpretation that fluid 
particles that are sitting just inside or outside the boundary surface when the star is unperturbed, will do 
so also when the star is oscillating. Hence no particles jump between the two phases, which in turn means 
that one is assuming that the transition reaction rate is much slower than the timescale associated with the 
perturbation. Since J- — implies continuity of ^, the Israel junction conditions are in this case satisfied 
simply by taking both ^ and Ap continuous. 

These two alternatives may be called fast and slow phase transitions, respectively jl7) . It is worth noting 
that the slow phase transitions are straightforwardly generalised to any type of matter; setting A_F = 
implies J-^ = which means that the junction conditions are satisfied by requiring ^ and Apr to be continuous. 
However, there is no obvious way to generalise the fast phase transitions to anisotropic matter. 

At the surface of the star where we wish to match our oscillating stellar model to the vacuum Schwarzschild 
solution with the same mass M as for the background solution, we set AF = since the material particles 
at the surface are of course staying on the surface during the oscillation. Thus F ~ and it follows that the 
physically relevant boundary condition is 



F 



Ap 



r{r — 2m) Ap 



(78) 



p' jr' 



[p + p) (m + ^nr'^p) 



Apr = 0. 



(79) 
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It also follows that <^ should be continuous, but since the frame vectors u°- and r° that are used to define A 
are not uniquely defined in the exterior vacuum region, this implies no further restriction at all. To make 
sure that we really get an oscillating solution of the same mass as the static solution, we should check that 
Am = at the surface. This will automatically be the case since 

Am = -iKrVC, (80) 

which vanishes when pr = Q. 

To our knowledge the general matching conditions (|74|l - (|76|l are new, and also the two above specialisa- 
tions to perfect fiuids seems not to be well known. 

At the centre we demand that the perturbed spacetime should be regular, which implies the vanishing 




leading to the condition 

hm = 0. (82) 



4 Elastic matter 

As described in Paper I we shall assume that the elastic matter equation of state is set up using a fixed 
metric as the sole structure on the three-dimensional material space X. Generically, this metric has to be 
spherically symmetric to be compatible with a spherically symmetric spacetime M , whether dynamic or 
static. Thus we may write this metric as 

kAB = tatb + tAB, (83) 

where Iab is the constant curvature two- metric on the SO (3) orbits and fA is the outwards directed unit 
normal to those orbits. Furthermore, compatibility with M being spherically symmetric also requires that 
the material space mapping ip should give an isometry between the unit curvature metrics on the SO (3) 
orbits of M and X, which in turn implies that the puUback tab — ip*tAB is 

iab = {f/rftab, (84) 

where f is the curvature radius of Iab- Since kab = 'ip*kAB should be orthogonal to the same clearly must 
be true for fa — 4'*^A, viz u°'fa = and fa — {r^fb)i'a- Ignoring the unphysical case when f is constant we 
may write fa — e^T>af for some function A of f. Thus it follows that f — and 

kab = {e^ffrarb + {f/rftab, (85) 
from which we read off that the linear particle densities (the square roots of the eigenvalues of A:";,) are given 

by 

Ur^e^f (86) 
nt — f/r. (87) 

Since nt is doubly degenerate we are, just as in the static case treated in Paper I, effectively dealing with 
a two-parameter equation of state p = Pcs{nr, nt) which comes from some full three-parameter equation of 
state p = p(ni, 77.2,713) by setting ni = n^, n2 = n^ = nt. The goal is here to specify the final perturbation 
equations ifSUIl and that are valid for a static background but with general matter to the case of elastic 
matter. Since the equations are expressed in terms of the matter variables p, Pr and g, it is convenient to 
work temporarily with n = nrn^ and nt as the two independent equation of state variables. The equation 
of state then gives Pr and q according to 

dp 

Pr = n- p 

on 



,.i„,|L. (89) 
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This in turn implies that the perturbation of p can be expressed as 



It foUows from eq. H87(l that 



Art Ant 

Ap= {p + pr) h 6q . 

n TLt 



Ant 



nt 



(90) 



(91) 



Substituting this form of Ant/nt into eq. (|9()(l and comparing the resulting expression for Ap with that 
implied by eq. (|50|1 . we can solve for An/n with the result 



An 
n 



1 

77 



2r' 



At this point it is convenient to replace ^ with the variable 

The expression for An/n then compactifies to 

An 



-2 1/ C 

r e — . 



This, in turn, implies that 



dpr dpr _2 V 

Apr = — — An + -— Ant ^ ~r e 
on ont 



dq 



dq 



Aq = — An+ —— Ant = -r n— — +7it- 

(jfl /'In . \ rirt -r*' rirt . 'r* 



dn 



Pt-: -6 q-n— - 



dq C 



C 



dn r' 



' dnt r 



where we have used the Maxwell type relation 

dpr 



dnt 



dq 



nt^r- = 6 n- q] , 



dn 



as well as the fact that the radial compressibility modulus Pr in the present variables is given by 

dn 



(92) 
(93) 
(94) 

(95) 
(96) 

(97) 

(98) 



We are now in the position to obtain a closed wave equation for ( by substituting eqs. H95|l and H96(l into eq. 
(|52(l . To display the result in a useful and compact way it is convenient to introduce the auxiliary variable 



where 



P 



Pre 



Y = 



Q - n— 
dn 



(99) 

(100) 
(101) 



Upon substituting eqs. H95|l and H96|l into eq. H52(l we now find that the radial perturbations of SSS elastic 
models are governed by the single wave equation 



(102) 
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where rj, P and Y are given above and 



W 



P 

V ^1 

r\\ 



Q = Qi + Q2 
P 



-21/ 



r\\ 





4r'\ 




v' [v' - 













p 



q[2v' - rit-P- — 

Ofit r 



(103) 
(104) 
(105) 

(106) 



with v^^^ = Pr/{p + Pr) being the squared speed of radially directed longitudinal (sound) waves. While it 
was convenient to use the variables (n, nj) to derive eq. H102(l . we shall now forget about them and change 
back to the variables (nr,nt) which makes Y and Q2 take the slightly different forms 



Y = 



6 



Q2^P 



dq 

drir 



3r' 



q ( 2v' ) + ( 2nr-z fit-^ — I — 



dnt 



Y' 



(107) 
(108) 



We also give the expressions for these quantities when the variables (n, z) are used, where z is the shear 
variable z = rir/nt- 



Y = 



6 



r(3, 



dq dq 

1 ^ " T ^TT 

on oz 



6 



q[2u'- 



3r' 



3z 



dq r' 
dz r 



Y' 



(109) 
(110) 



It is now easy to calculate Y and Q2 when using the particular quasi-Hookean equation of state used in 
paper I. The result can be displayed as 



Q2^P 



3r' 



q[2v' 2(/i + 3cr)- 



Y' 



(111) 
(112) 



where fx is the shear modulus, a the shearing energy density and = d(ln/i)/d(lnn), as described in paper 
I. Let us now abandon the covariant formalism and specialise to Schwarzschild coordinates (t, r) for the 
background. This is easily done, since we shall only need that the dot and prime derivatives of any scalar S 
are given by 

.dS „, ,dS 



S — e , S"— r , 
dt dr 



with 



We also separate out the time dependence according to 



c 



2m 



C(r), ^ e±''^*7?(r) 



which leads to the ordinary differential equation 

'^ + Y'n + {Q + uj''w)C = Q, 

dr 



(113) 
(114) 

(115) 
(116) 
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Clearly, this equation can easily be recast into the standard Sturm-Liouville form 



^(pf]+iQ + Lo'W)C = 0, Q = Q-PY'-^^. (117) 
ar \ dr J dr 

Comparing the definitions of the coefficient functions of the two equations (|116|l and (|117|) with those of 
the standard Sturm-Liouville equation that governs radial perturbations of SSS perfect fluids (c/. [121), we 
see that the functions P, W and the Qi part of Q are direct generalisations of the corresponding perfect 
fluid functions, while Y and the Q2 part of Q vanish identically when specialised to the perfect fluid case. 
Since eq. (|117|l is in Sturm-Liouville form whereas H116() is not, one could think that the former is the most 
useful way of writing the equation. However, we shall now argue that this is not the case. We saw in 
section 121 that ^ and Apr are contiuous variables if p and q are continuous for the background solution or, 
indeed, if we are considering slow phase transitions only. In that case the variables C and r] will by their 
definitions be continuous too. On the other hand, continuity of p and q does not imply that other matter 
quantities occuring in the coefficient functions of eqs. 1)116(1 and 1(117(1 are continuous. For instance, for the 
crust models considered in paper I both p and q are continuous all the way to the surface of the star but the 
shear modulus p. has a jump discontinuity at the interface between the interior perfect fluid phase and the 
exterior elastic one. This discontinuity in p propagates to other quantities such as the radial compressibility 
modulus Pr and hence the speed of sound v^w for radially propagating longitudinal waves. In particular, 
for these models the jump discontinuity of p in this case implies a jump discontinuity of PY , which in turn 
means that PdC,/dr = PY( + 77 should not be assumed continuous. Moreover 77 has a more direct physical 
interpretation than PdC,/ dr as it is directly related to the Lagrangian perturbation of the radial pressure 
and it is therefore ry and not Pd(/dr that should vanish at the surface of the star. For these reasons it 
is in general advantageous to work with eq. ((116(1 and to use C and r/, rather than ( and Pd(/dr which is 
common when studying perfect fluid starsjHI, as the two independent variables when recasting the second 
order differential equation pi6(l into a system of two first order ones. As a bonus, it will not be necessary to 
calculate the derivative d{PY)/dr that enters the expression for Q. The resulting first order system reads 

'^=Y(: + p-^ri (118) 
dr 

^--[(Q + c.2w-)C + Fr;]. (119) 
It should be supplemented by the boundary conditions at the centre = and surface r = R 

c 

lim - = (120) 
1 — >o r 

lim r] = 0, (121) 

r— »_R 

which follow from eqs. (jHU and ifT^ . 

We now wish to determine how ^ and 77 should behave near the centre. To do so we again employ the 
quasi-Hookean class of equations of state described in paper I. For this class it holds that 

P = Par-'^ + 0{r°) (122) 

Y = Yor-^ + 0(r) (123) 

W = Wor-^ + 0(r") (124) 

Q^Qor-^ + 0{r-^), (125) 
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where 



Po = e^''^{f3c + liic) >0, (126) 

^0 = o ^^4 - ' < yo < 3, (127) 

M^o = e''=(pc+Pc) >0, (128) 

go ^ -12e^-- <0. (129) 



The expansion coefficients Pq, Yq and Qo satisfy the relation 

Qo + PoYo{3-Yo) = 0, (130) 

which must be used to conclude from eq. 1)102(1 that ^ constant or C '--^ r"^ near the centre. The regularity 
condition ((82|l tells us that the latter behaviour is the physical one. We hence conclude that the variables C 
and 7] leave the centre according to 

C = Cr^ + 0{r'^) (131) 
77 = CPo(3-ro) + 0(r2)..., (132) 

where C is an arbitrary constant. Although we have derived this result for a restricted class of equations of 
states, one should recall that the shear is zero at the centre of the star and thus small in the vicinity of the 
centre. Hence the result should generalise to all equations of state with a physically reasonable behaviour in 
the limit of small shear, since the given class gives the most general physically reasonable behaviour in that 
limit. 

Now, the first order system consisting of eqs. 1(118(1 and ((119(1 is in general not a standard Sturm-Liouville 
system due to the presence of the function Y. Moreover the function Q has a non-standard r~* behaviour 
near the centre, which however goes over into the standard behaviour precisely when the leading order 
r^^ term of Y vanishes. Still, it is not difficult to show that the main Sturm-Liouville results for the standard 
case will remain valid for our somewhat more general setting, assuming continuity of C and rj. To see that all 
eigenvalues uj^ are real, let (CmVn) and {Cn',Vn') correspond to the eigenfrequences a;„ and ujn', respectively. 
Using eqs. ((118(1 and ((119(1 one readily finds that 

^ (Cnr/„' - Cw fin) = [iO^ - iO^, ) WUn' • (133) 

Integrating this equation over < r < R gives 

iCnVn' - Cn'Vn)\ = {Co^ - U:^,) / W'CnCn'dr (134) 

Jo 

Since the boundary conditions always make the function Cn'7Ti' ~ Cn'^jn (the generalised Wronskian) vanish 
at r = and r = i?, it follows from setting n' = n that — uj^ for all n, since W > for < r < i?. Since 
the eigenvalues thus indeed are real we can also choose the eigenfunctions Cn to be real and conclude that 

R 

WCnCn'dr^O IoxujI^ujI., (135) 



i.e. eigenfunctions corresponding to distinct eigenvalues are orthogonal with respect to the interval < r < i? 
and the weight function W . To sec that the eigenvalues are nondegenerate, form an infinite discrete sequence 
which is bounded from below but not from above, and can be enumerated such that 

ujI < < w| < . . . , (136) 
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with each corresponding C„ having exactly n internal {0 < r < R) separated nodes, it is convenient to make 
the Priifer type substitution 

C = Asinip, ri = Acos(p. (137) 
This leads to a first order differential equation for the phase ip, reading 

^ = p-^ cosV + Y sin 2ip+{Q + u?W) sinV- (138) 
dr 

Since tan(^ = QjT] ^ as r ^ 0, we can declare that should leave the centre according to 

p,./ ^y ^O{r% (139) 

which directly follows from eqs. H131|) and (|132|l . Using the observation that P is strictly positive on 
< r < i?, it is immediate that 

^ - P-I > 0, (140) 
dr 

at internal radii where is an integer multiple of tt, i.e. where ^ has internal nodes. This implies that the 
C = line of the 77-C plane is always crossed in a counterclockwise direction and that the nodes of C are 
separated since the origin can never be encountered unless we are dealing with the trivial solution ^ = 77 = 0. 
Based on the fact that W is also strictly positive on < r < i?, it is fairly easy to show that the phase, at 
each fixed internal radius r, is a continuous, monotonically increasing function of such that 

lim (/^(r, w^) = 00 (141) 

lim v'(r,w2) = 0. (142) 

— ^—00 

This is exactly what is needed since it implies that the phase function corresponding to C„ will then 
satisfy 

</'«(?'«') = "-'t^; n' = l,2,...,n (143) 

Jim <^„(r) = (n+ i) TT (144) 

where r„' is the radius of the n':th node of C„. The limit (|144|) corresponds to the boundary condition at 
the surface of the star being 

lim ry = 0. (145) 

) — 

Proofs of these statements on how (^(r, cj^) depends on uP' are presented in the appendix. The proofs are 
based on continuity for C and 77, but it should be remarked that they do not fail to be valid when the 
coefficient functions P, Y , Q and W have a finite number of jump discontinuities. 

The above results are very pleasing as they tell us that the radial perturbations of elastic stars is in 
essence a simple Sturm-Liouville problem that is a direct generalisation of the corresponding problem for 
perfect fluid stars. 

To conclude this section it is interesting to note that, like the perfect fluid case^HI, the eigenvalues may 
be derived from a variational principle. To this end we observe that the linear system H118(l - H119|) follows 
from the one-dimensional "time" -dependent Hamiltonian 

H = \[p-^Tf + 2YCri+{Q + oJ^W)C^]. (146) 

Indeed, since 



dC _ dH 

dr dr] 
drj dH 
dr dC, ' 



(147) 
(148) 
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we see that 77 can be interpreted as the conjugate momentum of C. From the corresponding Lagrangian, 



under variations of C that keep the boundary conditions satisfied, are precisely the eigenvalues uj^. 
4.1 Stability and the mass-radius curve 

A technique which has proven very useful for stability considerations of perfect fluid SSS stars is based on 
the analysis of static configurations of a whole family of stars corresponding to a specified equation of state. 
Given appropriate boundary conditions, for a given central pressure pc, there exist a unique solution to the 
TOV equations, with a corresponding gravitational mass M(pc) and total (Schwarzschild) radius R{pc)- The 
set of equilibrium configurations of a particular EOS can thus be represented as a curve in the M-R plane. 
If shells of infinite density are disallowed the curve will be continuous. Moreover, in the absence of first 
order phase transitions the curve is guaranteed to be C^ISOI (viewed as a submanifold of the M-R plane). 
Postponing for the moment the discussion of possible implications of different kinds of phase transitions we 
indicate how the M-R curve can be used to determine the stability of a given star. The original arguments 
were found by J. A. Wheeler in the context of cold white dwarf and neutron star models and are presented 
in IE] (but see also PT]V 

As a first observation, note that since the radially perturbed models are matched to the same exterior 
Schwarzschild solution as the corresponding background models, it follows that static perturbations {lo^^ = 
for some n) are allowed only at extremum points [dM/dR = 0) of the M-R curve, since only at such 
points can static models be perturbed to static neighbours of the same mass. Clearly then, owing to the 
continuity of the eigenfrequencies viewed as functions of the central pressure, the stability of the star can 
only change at such extrema. Moreover, general arguments following from the Sturm-Liouville properties of 
the perturbation equation show that the mode changing stability at the extremum is even (n = 0, 2, . . .) if 
dR/dpc < and odd {n = 1, 3, . . .) if dR/dpc > 0. For equations of state whose low pressure behaviour is 
known to give stable configurations (this obviously includes neutron stars since small iron balls are stable), 
it follows as a corollary that a mode is (de)stabilised if the curve turns (counter)clockwise at an extremum. 

Due to the general nature of the arguments they can be directly applied also to the elastic case once 
the uniqueness and Sturm-Liouville properties have been established. It should be clear from the results 
of Paper I (see also 22, 23 ) that, in the absence of jump discontinuities in the energy density and/or the 
tangential pressure, the equations of stellar structure possess a unique solution for each central pressure. 
Although no proof is presented here it is fairly obvious that the M-R curve is under these restrictions. 
Moreover, the Sturm-Liouville properties have been rigorously demonstrated in this case so we conclude that 
the static stability analysis is straightforwardly generalised to the elastic case. 

However, in situations with first order phase transitions in the density and/or the tangential pressure, 
the situation is more subtle. For perfect fluids only the energy density can be be discontinuous implying that 
the matching is uniquely determined by the junction conditions. However, when there is a transition to an 
exterior elastic phase, the tangential pressure just outside the phase boundary is not specifled by the junction 
conditions, as was shown explicitly in paper I for the case of a perfect fluid interior to a quasi-Hookean type 
elastic phase. This arbitrariness can be removed by considering the value of the pressure anisotropy to be 
given by a parameter specified as a part of the equation of state. If this view is taken, each central pressure 
corresponds to a unique solution. In the presence of such first order phase transitions the M-R curve is well 




(149) 



and standard methods |19|. it can be readily shown that the extremal values of the quantity 




(150) 
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defined, continuous, but may not be C^. For the crust models presented in Paper I with M-R curve and 
compactness plots, the arbitrariness was removed by simply setting the pressure anisotropy to zero on the 
outer side of the phase boundary, which resulted in models with continuous energy density and tangential 
pressure. Thus the static M-R curve approach to stability can definitely be used for these models. 

Had we instead made other choices for the pressure anisotropy (in fact we did play around with other 
choices, as mentioned in the paper) the situation would be more unclear since 1) the M-R curve may not be 

and 2) we have not proved that we are dealing with a Sturm-Liouville type problem. The same is true 
for the models with elastic cores that were also presented in Paper I, as well as for more general models for 
which p and pt are not continuous. 

For a perfect fluid with a jump discontinuity in the density at some transition pressure pi, the behaviour 
of the M-R curve was studied by Lindblom|2U|. Parametrising the discontinuity by A = (p+ — p-)/{p- 
where p± are the densities just above and below the phase transition, he found that there exists a critical 
value of this parameter Ac = (p- + 3pi)/2(jO_ +pi) where the slope of the M-R curve has a discontinuity. 
Beyond the critical point the curve actually turns 180°, leaving its tangent continuous. A similar analysis 
as Lindblom's for elastic matter should be very useful, but is postponed for later work. 

Of course, if one is uncertain about the stability of a given model one should always calculate the 
eigenvalues. For the relativistic polytropic stellar models with quasi-Hookean elastic crusts considered in 
Paper I we computed the first few modes and confirmed that the fundamental mode did indeed have a zero 
at the maximum of the M-R curve. For the core models, we do not know theoretically that instability 
should set in at the maxima of the M-R curves, since the tangential pressure and energy density are in 
fact discontinuous at the boundary. However, assuming that the phase transition is slow, we find that these 
models become unstable at, or at least very close to, the mass maximum. 

5 Conclusions / discussion 

With the results of the present paper at hand we are now able to numerically determine the frequencies 
of adiabatic radial perturbations of general relativistic SSS stars with elastic matter sources, without much 
further ado than in the much studied special case of perfect fluid stars. For instance, both numerical 
schemes employed by Kokkotas and RuofflH] (based on the shooting method and the Numerov method) 
are very easily modified to allow for elasticity. Moreover, we established that the mass-radius curve, under 
certain conditions, can be used to assess stability without having to calculate a single frequency. 

As expected, the spectrums of the radial perturbations of elastic SSS models with moderate shear moduli 
do not differ much from the corresponding spectrum when the shear modulus is identically vanishing, as is 
the case for perfect fluid models. 

Although a thorough treatment of general phase transitions is postponed to later work, it seems clear that 
modest transitions do not have a large impact on the stability properties. The connection between large 
phase transitions and stability has, to our knowledge, not been conclusively established even for perfect 
fluids. It would be very useful to determine the precise relation between the behaviour of the M-R curve 
and the sign of the eigenfrequencies, for slow as well as for fast phase transitions. We feel that it would 
probably be wise to treat the perfect fluid case first, since the elastic case is complicated by the additional 
freedom of specifying how the transition is to be made. It deserves to be mentioned that, as far as the crust 
of a real neutron star is concerned, the density discontinuities are not thought to have A larger than about 
10~^Ac|211. However, exotic phases in the cores of neutron stars might lead to larger discontinuities. An 
extreme example is provided by a strange quark star with a normal matter envelope, in which case A can 
be as high as 10^ Ac or more[25j. 

Since we now have the spherically symmetric equilibrium configurations as well as their stability under 
control, we are ready to take on more complicated, and hence more interesting problems. The first step 
is of course to treat non-radial first order perturbations of various kinds. Work on time-dependent axial 
perturbations of prestressed configurations is already on the wav|26| applying the general formalism of 
Karlovinij27|. We also plan to perturb SSS models into stationary neighbours, by setting them into rotation 
as well as changing their shapes, the latter being impossible for perfect fluid stars (to first order). What 
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remains as far as first order perturbations arc concerned is the time-dependent polar sector. However, unlike 
the situation in the axial sector, the matter and metric perturbations (to the extent that they can be talked 
of separately) are coupled in a highly non-trivial way already for perfect fluids, which would probably make 
the effects of elasticity very hard to assess. It therefore seems much more feasible as well as interesting to 
instead directly treat rotation to second order in the spirit of Hartle and Thorne f^ I20| . Indeed, such a 
setting would provide a sound basis for discussing very interesting astrophysical phenomena such as free 
precession, quakes and glitches of neutron stars, all of which are observed in pulsars PU]. 

Although of less astrophysical interest, it would be fun to find some exact SSS solution with a physically 
reasonable equation of state that is found to be stable using the results of this paper. We have in fact 
recently found a few solutions with regular centres, finite radii, everywhere causal wave propagation and 
positive energy densities and eigenpressures|3J|. One of these solutions has the central pressure as a free 
parameter and does consequently possess a M-R curve. Using the formalism of the present paper we found 
that instability sets in at the maximum of the AI-R curve as expected. 

Appendix 

In this appendix we prove that the phase function ip, satisfying the differential equation 

'^ = p-^cos^tp + Ysm2tp+{Q + \W)sm'^ip^:F{r,tp,X), X^uj^, (151) 
dr 

is a continuous, monotonically increasing function of A that satisfies 

lim <y9(r. A) = oo (152) 

A— *oo 

lim ip{r,X) = 0. (153) 

A — * — oo 

Since eq. (|151|l is singular at r = and possibly also at r ~ R, we must restrict ourselves to prove the 
statements for values of r in the open interval < r < i? only. 

The continuity directly follows from the continuity of F{r, A). To prove monotonicity we observe that 
the second term in the expansion H139() for Lp is 

(154) 

Thus it follows that given Ai and A2 with Ai < A2, the difference between the corresponding phase functions 
Lpi and will behave according to 

near the centre. Let > be a radius sufficiently small for fi{ri) < ip2{f'i) to hold. This gives us the two 
first order differential equations 

^^ = nr,^i,Ai) (156) 

^ = nr,^2,A2) (157) 

(158) 



subject to the initial condition inequality 



vi{n) < Mn). (159) 



Since 

i^(r,<P2,A2) > F(r,v?2,Ai) forO<r<i?, (160) 
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it follows from a standard theorem on differential inequalities that 

ipi < ip2 for 0<r<R, (161) 

given that F{r, ip, A) is Lipschitz continuous in ip and continuous in r on this interval. The former condition 
is clearly satisfied since F{r, Lp, A) has a simple trigonometric dependence on Lp. The possible complication 
that F{r, ip, A) may have a finite number of jump discontinuities in r of the type 

F_ = lim_F(r,v5,A) ^ F+ lim F(r,v3,A), |i^+| < oo, < n < R, (162) 

r — >r. 1 — 

fa h 

is not a problem. Indeed, since we assume that ip is always to be matched continuously across such radii 
it is very easy to show that the theorem extends to a situation when F{r, (j), A) has discontinuities in r of 
this type. We have thus proved that ip{r, A) is a strictly increasing function of A for all r e (0, R). 

To prove the limit (|152|l we begin by claiming that it is not difficult to convince oneself that given the 
expansions H122|l - H125|l subject to the relation H13U|) . there exist constants > and A^ > for which the 
following statement is true: If and Wq are any constants satisfying 

< P^^ < Pq-i (163) 
< 1^0^ < Wq (164) 

then for all r G = (0, r^], for all ip and for all A > A^, 

F{r,ip,X) > P-^-^ cos^ip + sm2ip + {Q'' + XW^) sin^ip =: F^(r,(/),A), (165) 

where 

= Yor-^ (167) 

= W^r-^ (168) 

Q^ = Qy^ Qi = -Po^Y„{3-Yo). (169) 
This will be true if we can choose and A-'- such that, for all r <E I^, 

p-^-P^-^>0 (170) 

Q-Q^ +X{W ~W^)>0 (171) 

(p-i _ pi-i)[g - + X{W - W^)] -{Y- Y^f > 0. (172) 

To achieve this we first note that 

P-^ - P^-i = (Po-i - Pt^y^ + Oir^) > (173) 
W -W^ ^ {Wq -Wi^y-^ +0{r°) > (174) 

can be made valid on by choosing sufficiently small. Furthermore, since the coefficient of 

Q-Q^ = {P^- P)Yo{3 - Yn)r-^ + 0{r-^) (175) 

is positive (or zero if Yq = 0), we can choose A^ sufficiently large for the inequality (|171|l to hold. It now 
follows directly from Y — Y^ ^ 0{r) that the remaining inequality (|172|l can also be made valid, since it 
shows that it suffices to set A^ to an even larger number, if necessary. Let us now consider the two differential 
equations 

^ = F(r,^,A) (176) 
ar 

^^FHr,^K^), (177) 
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and focus on the unique solutions that leave the centre according to 

V^JTTT^y + Oir'') (178) 

Since 

we can always, for every fixed value of A, find G such that 

^{n,\)> ^^{n,X). (181) 

It then follows from eq. ptiSf) and the already used theorem on differential inequalities that 

^{r,\)>^^{r,\), (182) 
for all r ^ and all A > A^. We now use that (p^{r, A) can be found in closed form and is given by 



1 sin(fcr) - A:rcos(fcr) n , W} \ , „„, 

(kry sm(fcr) - Y^)[sm[kr) - kr cos{kr)\ y p^^^ 

Since we can hence infer that 

lim ip^{r, A) = oo (184) 

A — ^oo 

for every r > 0, the same clearly is true for (^(r, A) for every r g Finally, we only need to note that since 
(/3(r, A) can never cross a line (f = tt-tt downwards, it follows that the result extends to the whole interval 
< r < i?. 

Now, to prove the limit H153|l we start off by instead bounding F(r, </), A) from above near the centre, 
much in the same way as when we bounded it from below by i^^(r, 0, A). We thus take two constants Pq 
and Wq such that 

Pi' > P^ > (185) 
> M^o > 0. (186) 

Then there exist > and A^ < such that for all r G = (0,m^], for all ip and all A < — |AT|, it is true 
that 

F(r,(^,A) < p1^"icosV + 5^^sin2(p+ (Q'f + AVKT)sinV =: F^r,(f),X), (187) 

where 

= (188) 

= Yor-'^ (189) 

H/T ^ w^-'^ (190) 

W^=Qlr-^ Ql = -P^Yo{3~Yo). (191) 

Letting Lp{r,X) and 93^ (r. A) be the unique solutions of 

^=P(r,0,A) (192) 
dr 

^=FHr,<l>,X), (193) 
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that leave the centre as 

V{r,\)^ ^ r^ + 0{r^) (194) 
fQ{6 — Joj 

vHrA)= / ^^ r^ + 0{r^), (195) 
(-^ - ^o) 

then, if r e and A < -|AT|, 

ip{r,\) <if>\r,\). (196) 
For negative A the exact form for {r, A) is given by 



t, 1 sinh(A:r) — /cr cosh(fcr) o , /w^n|A| , , 

^ ' ' pj (fcr)2 sinh(fcr) + ro[sinh(fcr) - fcrcosh(A:r)] V 

from which we can see that 

Hm tp\r, A) = for all r > 0. (198) 

A — ' — oo 

The same is hence true for (p(r, A) for every r £ /^^ since it cannot cross the line = downwards. Let us 
now choose A sufficiently negatively large so that (l){r'^) < e for some e £ (0,7r). Let us also assume that 
(p{rQ,X) = e for some tq G {r\R). However, since sine 7^ we can make dip/dr negative at r = tq by 
choosing |A| even larger, if necessary, implying that ^p{r, A) cannot reach the line (p — e from below if A is 
sufficiently negatively large. Thus we have proved the limit H153|l since e can be chosen arbitrarily small. 
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